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We present a detailed model study of exciton transfer processes in donor-bridge-acceptor (DBA) 
systems. Using a model which includes the intermolecular Coulomb interaction and the coupling 
to a dissipative environment we calculate the phase diagram, the absorption spectrum as well as 
dynamic equilibrium properties with the numerical renormalization group. This method is non- 
perturbative and therefore allows to cover the full parameter space, especially the case when the 
intermolecular Coulomb interaction is of the same order as the coupling to the environment and 
perturbation theory cannot be applied. For DBA systems up to six sites we found a transition to 
the localized phase (self-trapping) depending on the coupling to the dissipative environment. We 
discuss various criteria which favour delocalized exciton transfer. 

PACS numbers: 71.27.+a Strongly correlated electron systems, 71.35.Aa Frenkcl Excitons and Self- Trapped 
Excitons, 71.35.Cc Intrinsic Properties of Excitons; Optical Absorption Spectra. 



INTRODUCTION 

Exciton transfer belongs to the key processes in many 
chemical and biological systems, organic based nano- 
structures and semiconductors 0, 0, 0]- The progress 
in manufacturing molecular electronic devices, biological 
hybrid systems, and model systems based on quantum 
dots, nanoscale molecular aggregates and bio-engineered 
proteins opens the door to understand these fundamen- 
tal processes |4| and also to find applications in (bio-) 
molecular electronics, biosensing, and quantum compu- 
tation 

Excitons are electron-hole pairs which do not transfer 
charge but energy by deexciting a donor molecule fol- 
lowed by the excitation of an acceptor molecule. The 
radiationless excitation transfer is caused by dipole and 
exchange interactions and proceeds via a short lived vir- 
tual photon In this work we consider Frenkel ex- 
citons where the exciton is a molecular excitation 
with an electron in the lowest unoccupied molecular or- 
bital (LUMO) and a hole in the highest occupied molec- 
ular orbital (HOMO) on the same molecule. Here the 
Coulomb coupling of the electron-hole pair is much larger 
than the hopping matrix clement of a single hole or elec- 
tron. Exciton transfer where the Frenkel exciton con- 
cept can be applied, occurs in many bio- molecules, e.g. 
rhodopsin, porphyrins, blue copper protein, carotenoids, 
and chlorophylls. A well studied molecule is the light 
harvesting antenna (LH II) from the bacterial photosys- 
tem Rhodopseudomonas Acidophila. It is characterized 
by a symmetric structure and composed of nine identical 
units forming a ring. Each unit is composed of a chloro- 



*present adress: Department of Physics, Renmin University of 
China, Beijing 100872, China 



phyll dimer. The light-harvesting complexes store and 
transfer excitations with high efficiency. 

In the photosynthetic process a LH-II ring absorbs a 
photon. The excitation is transferred to other LH-II rings 
and sent via the LH-I ring to the reaction center and 
then converted to chemical energy. The excitation in the 
LH-II ring B850 can move over the whole ring - it is de- 
localized over the ring. In other rings, such as B800, the 
excitations are usually considered to be more localized 
4J. Furthermore, mechanisms exist which dissipate exci- 
tation energy to safe the organism from damage [s| . The 
degree of derealization depends strongly on the coupling 
to the vibronic environment and may be crucial for the 
function of the specific protein. 

The interpretation of optical spectra J9( requires a the- 
ory which incorporates both static and dynamic disorder. 
If the fluctuations of the protein environment occur on a 
much larger time scale than those of the excitonic system, 
the disorder is regarded as static. Such a static disorder 
can be treated by a thermodynamic average. The dy- 
namic disorder stems from the coupling of the electronic 
degrees of freedom to the fluctuations of the environment. 
In the present paper we will study the effect of dissipation 
while neglecting static disorder. 

A full ab initio quantum chemical calculation of 
molecules which show exciton transfer reactions is impos- 
sible; therefore it is reasonable to investigate the system 
using simple models which, nevertheless, cover the rele- 
vant physics of the problem. The most elementary non- 
trivial model which describes quantum dissipation is the 
well studied spin-boson model 0, . It can be viewed 
as an archetype for modelling the system-environment in- 
teraction in bio-molecules in which the electronic degrees 
of freedom couple to a dissipative environment. 

A variety of theoretical methods have been devel- 
oped to calculate absorption spectra and rates (see e.g., 
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EHQHEGllEH)- Some investigations of ex- 
citon transfer systems were based on perturbation theory 
in the exchange coupling between the excitons or in the 
coupling of the electronic system to the vibrations. If the 
exciton-vibrational coupling a is weak compared to the 
dipole-dipole coupling density-matrix theory is used 
If the intcrmolccular Coulomb coupling J is small we are 
in the limit of nonadiabatic exciton transfer. Here per- 
turbation theory is applicable which leads to the Forster 
equations Q. 

A key challenge for a theoretical study of exciton trans- 
fer is to cover the whole range of possible behaviour, from 
coherent to incoherent transfer or even localization or 
self-trapping of the excitations. Here we use the non- 
perturbative numerical renormalization group (NRG) to 
calculate equilibrium properties of the exciton system in 
the full parameter space. We give a detailed study of 
the phase diagram, dynamic equilibrium properties for 
chains and rings up to six sites, and the frequency de- 
pendent linear absorption spectra of excitons in a dimer 
and trimer molecule as a function of the coupling to the 
bosonic bath. The behaviour of the system is governed 
by the competition between the couplings J and a which 
determine whether the excitations are delocalized or lo- 
calized. 

The models we are considering here describe general 
electron-boson systems with a limited number of quan- 
tum states on a few sites with a coupling to a (quan- 
tum) dissipative environment. We restrict ourselves to 
donor, acceptor, and bridge molecules with only two elec- 
tronic levels per molecule, neglecting the spin degree of 
freedom. If the flux of photons is sufficiently low then 
the exciton migration in systems such as a pigment net- 
work can be satisfactorily modelled by a single excitation. 
We show that in the single-exciton subspace the multi- 
site electron-boson model can be mapped to a multi-site 
exciton-boson model. The two-site exciton-boson model 
is identical to the spin-boson model. The exciton system 
is coupled to all degrees of freedom of the (protein) envi- 
ronment which is modelled by an infinite set of harmonic 
oscillators. After the discussion of the various models 
in Sec. II, we introduce the NRG approach used here 
in Sec. III. Section IV is devoted to the results for the 
phase diagram, dynamic properties, and the absorption 
spectrum. We show how the degree of derealization de- 
pends on the different Coulomb interactions, the coupling 
to the bosonic bath, and the geometric structure. 



MODEL 

In general we describe the problem by a small elec- 
tronic system like a short chain or ring of molecules with 
the electronic part H e \ coupled via the -f/ci-bath part to 



the vibronic degrees of freedom incorporated in -ffbath^ 

H = H c \ + iid_bath + -ffbath ■ (1) 

The simplest possible model for the electronic part is 
a two-state system with the two states corresponding to 
the electron being located at the donor or at the acceptor 
site. In this case, the electronic part can be modelled via 

H c \ = ^ £ i c i c * ~ K c a Cd + c^ca) ■ (2) 

i=A,D 

Apparently, models of this kind where the electrons are 
allowed to hop between donor and acceptor sites (with 
hopping matrix element t) are connected to electron 
transfer problems. In this paper we do not consider such 
hopping processes and focus on excitation transfer in- 
duced by a two-particle interaction term in the Hamilto- 
nian. The electronic part of the Hamiltonian then takes 
the following general form 

H el = J2H(i,l,<r)+ V(i,j;k,l;a,a') , (3) 

where the first term represents the local part with on-site 
energies, on-site Coulomb interactions as well as possible 
spin-orbit couplings. The second part describes all non- 
local terms, mainly the non-local Coulomb interaction 
(as a possible extension also the single-particle hopping 
term). In principle, the parameters of the Hamiltonian 
eq. @ can be extracted from quantum chemical calcu- 
lations [iolEH. 

In the current work, we explore excitons in a chain or 
ring using the Hamiltonian 

H ci = ^ e ik4k°ik + J ^4k c ]l c 3l' c ^' ■ (4) 

ik ijlkl'k' 

The operators CyJ denote annihilation (creation) opera- 
tors for the electrons on site i in the level k, Cik is the on- 
site energy and Jij the exchange interaction between site 
i and j . The electronic part of each level is represented by 
a two-level system. We neglect the single-electron hop- 
ping between neighbouring sites as well as spin degrees 
of freedom. Therefore, the excitations (Frenkel excitons) 
can only be transferred via the Coulomb coupling. Sim- 
ilar models were discussed, e.g., in [a. IwL llSl 12^ 

The Hamiltonian eq. I@J describes a coherent motion 
of the excitation through the whole system; this coher- 
ence can be destroyed in the presence of a dissipative 
environment. Here, the coupling to the environment is 
due to the change of the dipole moment of the molecule 
during the transition. Simulations showed that the 
coupling involves essentially all nuclear degrees of free- 
dom of the protein which have to be described quantum 
mechanically. Even at physiological temperatures there 
are many degrees of freedom in proteins with frequencies 
high enough to make a quantum mechanical description 
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necessary |23J . We represent the vibrations of the envi- 
ronment by a set of harmonic oscillators similar to the 
spin-boson model. 

The last term in eq. describes the free bosonic bath 



J DA 



1 >a 111 



(5) 



with the bosonic annihilation (creation) operators bn^ ■ 
The second term in eq. £Q| describes the coupling of the 
electrons to the bosonic bath 



1-bath 



gin™ Y A ™ ( b i + K 



(6) 



where A„ is the coupling strength to the nth oscilla- 
tor. We consider here a dipole coupling so that the sum 
Y^i 9ifii is the polarization operator of the electronic sys- 
tem. The values of the parameters gi will be specified 
below. 

In analogy to the spin-boson model 0, [n| , the cou- 
pling of the electrons to the bath degrees of freedom is 
completely specified by the bath spectral function 



(7) 



Several parametrizations of J{u>) have been studied in 
the literature 01 • For a given system, the bath spectral 
function can also be calculated using molecular dynamics 
simulations 0|. Here we restrict ourselves to a simple 
ohmic spectral function and will use more realistic spec- 
tral functions in a future study. 



Dimer 

For the dimer (that is the two-site electron-boson 
model as sketched in Fig. ^) the Hamiltonian reduces 
to 

#dimcr = H e \ + -f/el-bath + -Hbath , (8) 

with the electronic part 

H c i= Y £i4°i + J ( c £» c Ai c ^i c -Aq + /l - c -) '( 9 ) 
i=D ,Di,A ,Ai 

and the coupling term 

-^el-bath = 



Y 9imJ2 X "( b n + b n) ■ (10) 

i=Do,Di,Ao,Ai n 

The indices Dq,Aq indicate the ground state on the 
donor/acceptor and D\,A\ the first excited states on 
the donor/acceptor. For the gi we put go t = —gD = 
~9 Ax = 9a = \- 

The electronic degrees of freedom in this subspace can 
be represented by the four-dimensional basis 
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FIG. 1: Schematic view of the two-site electron-boson model. 
The transfer integral for the exciton transfer between donor 
(D) and acceptor (A) sites is given by J. Dissipation in the 
exciton transfer process is due to the coupling to a common 
bosonic bath. 



with the notation |D;A) describing the donor/acceptor 
in the ground state (D/A = 0) or in the excited state 
(D/A = 1). Introducing the notation 



(12) 



we arrive at the matrix M = Ma + Mb = (i|i?dimer|i) 
= 1..4), where the matrix elements are taken only 
with respect to the electronic states: 



M 



and 



/ e D „ + CA 

e Dl + e Ao 
J 





J 

e_D + e Al 

e Dl + e Al J 



M h 



(Y \ 

Y-X 

Y+X 

V Y J 



With e Do — e Ao — and e Dl = e D , e Al 
eigenvalues of the electronic part M are 



(13) 



(14) 



e A the 



E 1 = 
Ei — 

^3,4 

■ 



, 

e_D + (-A , 



(e D - e A ) 



J 2 



|i) = {|0,0),|l,0) ) |0,l),|l,l)} , 



(11) 



The eigenstates with energies £3,4 are linear combina- 
tions of the basis states |1, 0) and |0, 1). 

The Hamiltonian of the dimer eq. (jHJ) can be decom- 
posed into subspaces of zero, one, and two excitons. For 
the subspace with one exciton in the dimer, the basis 
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FIG. 2: Schematic view of the three-site electron-boson 
model. The Coulomb matrix element of excitons between 
donor (D), bridge (B) and acceptor (A) sites is given by J. 
Dissipation in the exciton transfer process is due to the cou- 
pling of the electronic degrees of freedom to a common bosonic 
bath. The excitation transfer shown in the figure is due to 
the coupling Job 

DBA _ Jd b DBA 
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FIG. 3: Schematic view of the three-site (chain) electron- 
boson model. 



eq. jnj reduces to \i) = {|1, 0), |0, 1)} and the Matrix M 
reads 



,1/ = 



e A + Y - X J 

J e D +Y + X 



(15) 



This matrix allows for an exact mapping onto the spin- 
boson model (for a similar discussion, see Ref. l24j) . The 
model is equivalent to eq. H22|) with N = 2. 



Trimer 



The Hamiltonian for a donor-bridge-acceptor system 
in a trimer geometry (see Figs. [21 and |3J) takes the form 



= E 

i=Z»o,Di,Ao,Ai,B ,Bi 



£j, C„- Ci 



+ Jda {c Dq c s Ai c Di c Aq + h.c^j 
+ Jab {c i Bo c\ i c Bl CA a + h.c 
+ Jbd [ c b c ' { d 1 c BiCDq + h.c 



(16) 



The electrons are coupled to the bosonic bath via: 
#el-bath = E 9ini E Y ( 6 " + 6 «) 

i=D ,D 1 ,A ,Ai,B ,Bi n 

(17) 

For the gi we choose gu l = —QDx = — 9Aq = 9Ai = 1 an d 
9B = 9b 1 = 0. 

The basis for the electronic degrees of freedom is now 
composed of eight states: 

|1) = |0,0,0),|2) = |1,1,1>,|3) = |1,0,0),|4) = |0,1,0), 
|5) = |0,0,1),|6) = |1,1,0),|7) = |0,1,1),|8) = |1,0,1), 

where the 0(1) indicates an occupied first (second) level, 
respectively, with the notation \D,B,A) for the occupa- 
tion of donor, bridge, and acceptor molecule. A direct 
hopping of the exciton is possible from the donor to the 
acceptor or to the bridge and from the acceptor to the 
bridge and back (see Fig. El- The matrix-elements now 
read: 



<ll#ll> 


- £d + es 


+ £A - 


h-F-1 




+ 5D - 


^9Ao)X, 


(2|H|2) 


= £Di + £Ai 


+ - 


hf-1 




+ 9D1 - 


>rgBl)X, 


(3|#|3) 


= e^o + esc, 


+ - 


hF4 




+ 9 B - 


\-9D t )X, 


(4|i2]4) 


= CD + e Bi 


+ £A ~ 


hF-1 




+ 9 D ~ 


^9A )X, 


(5\H\5) 


= CDo + e Ai 


+ e B ~ 


hr-i 


" (5D 


+ 9A t ~ 


\-9B )X, 


(6\H\6) 


= e,4i + £Bi 


+ e Dl - 


i-y-r 




+ 9A 1 - 


Vg Dl )X, 


(7\H\7) 


= £D + e Bi 


+ e Al - 


hf-i 


■ (ffBi 


+ 9 Do - 


hg Al )X, 


(8\H\8) 


= £Dj + €b 


+ e Al - 


hf-i 


- (5Bo 


+ 9D 1 - 


hg Al )X, 


(4\H\5) 


= (5|ff|4) = 


Jab, 










(6\H\8) 


= (8|fl-|6) = 


Jab, 










(5|i2]4) 


= (4|JT|5) = 


Jbd, 










(7\H\8) 


= (8|^|7) - 


Jbd, 










(3\H\5) 


= (5|i?|3) = 


Jad, 










(6\H\7) 


= (W) = 


Jad- 











For e A0 = eso = £_do = e and J AB 
the eigenvalues are given by 



Jbd 



(18) 
Jad = J 





= o, 




E 2 


= 3e, 




E 3 


= e — 


J, 


Ei 


= e + 


2J, 


E 5 


= 2e 


-J, 


E 6 


= 2e 


■f 2J. 



For the chain as in Fig. we set J^u = and J^s 
Jbd = J ■ The resulting eigenvalues are: 



0, 
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E 2 = 3e, 
E 3 = 2e, 

E 5 = e±V2J, 
E 6 = 2e±v / 2J. 

In the subspace with only one exciton, the basis consists 
of the three states |1, 0, 0), |0, 1, 0), and |0, 0, 1), and the 
matrix reduces to 

(e D + Y + 2X J AB J ad \ 

M = Jab e B + Y J BD . (19) 

V Jad Jab e A + Y-2X/ 

The eigenvalues for the trimer as in Fig. [21 are 

E-l = e-J, 
E 2 = e + 2J, 



and for the chain as in Fig. [3| 

Ex = e, 

E 2 = e±V2J. 

The model is equivalent to eq. 11' 21 with N = 3. 

Multi-site exciton-boson model 

In the single-exciton subspace, the fermionic degrees of 
freedom of the models introduced above can be mapped 
onto operators a\ for hard-core bosons corresponding 
to the creation and annihilation of an exciton at site i. 
This results in general multi-site exciton-boson models 
with N sites defined by 



-ffmulti — H K + -f/x-bath 

with the electronic part defined as 

N 



bathi 



H x — ^ Jija 



(20) 



(21) 



The parameters for i =/= j are the transfer integrals 
between site i and j. As before, we only consider nearest 
neighbour interactions. The diagonal elements J a are the 
on-site energies £j at site i. We perform a constant shift 
of the Hamiltonian by Ja = Jjj{^i,j) and arrive at 



N 



(22) 



For the coupling term we assume the following form 

AT 

# x -bath = 9ia\ai Y K (bl + K) (23) 



with gi = (i - (N + l)/2). For N = 2 the model is 
equivalent to the spin-boson model with the matrix M 
as in eq. (|15|l and for N — 3 equivalent to eq. 1)19(1 . 

The eigenvalues of H x for N = 4 are B = 0, ±2 J 
for the ring geometry and E = ±J for the chain. The 
eigenvalues for the ring with N = 5 are E = 2 J, (V5 - 
l)J/2, (— v/5 - l)J/2 and for the ring with iV = 6: 
S = ±J,±2J. 



METHOD 

The models we are considering here are completely 
specified by the parameters of the electronic system and 
the spectral function J(ui) (defined in eq. J7J) which can 
be estimated in a classical molecular dynamics simula- 
tion. We are using here an ohmic form: 



J (ui) = 2naujQ(u> — lu c ) 



(24) 



where a is the dimensionless coupling for which we use 
values in the range 0.01-2. The parameters J and t are 
all measured in units of lo c . Typical values of Hlu c are of 
the order of 1 to 10 meV. 

As described in the introduction, basically all degrees 
of freedom of the bosonic bath (the dissipative environ- 
ment) are relevant for the behaviour of the electronic or 
excitonic system. So it is not possible to disregard high 
energy states even if we are interested in low tempera- 
ture properties like the coherent behaviour for temper- 
atures smaller than the characteristic temperature T*. 
The renormalization group ansatz is designed for prob- 
lems where every energy scale contributes and perturba- 
tion theory typically shows logarithmic divergencies at 
small frequencies (energies) when the temperature goes 
to zero. 

In order to keep the paper self contained we explain 
the Numerical Renormalization Group (NRG) method 
for the bosonic bath in detail. 

Originally the NRG was invented by Wilson for a 
fermionic bath to solve the Kondo problem |2Jj|, |2g . The 
fermionic NRG is a standard and powerful tool to in- 
vestigate complex impurity problems with one or more 
fermionic baths. Only recently, the method was extended 
to treat quantum impurity systems with a coupling to a 
bosonic bath [2?lE8j |. Here we focus on equilibrium quan- 
tities (recently it was shown that the NRG can also be 
applied to non-equilibrium situations H^|.) 

For the numerical renormalization procedure we start 
from the Hamiltonian written in a continous form (see 
the detailed discussion in Ref. |26 



H = ff d + J deg(e)b% + fcftJ J de(b\ + b e ) 



(25) 
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FIG. 4: Logarithmic discretization of the spectral function 
J(w). 
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FIG. 5: Scheme of the bosonic chain-NRG. The boxes repre- 
sent the iterative diagonalization. 



The function g(e) is the dispersion of the phonon bath 
and h(e) characterizes the coupling between the elec- 
tronic system and the bath. Both functions are related to 
the spectral function J(uj), see Ref. l28i We start defining 
the renormalization group transformations by a logarith- 
mic discretization of the bath spectral function (Fig.^J in 
intervals [A~™ +1 ,A _n ], with n = 0, 1, . . . , oo and A > 1 
the NRG discretization parameter. The discretization is 
exact for A — ► 1 and works still very well for A = 3. (Here 
we are using A = 2) 

Within each of these intervals only one bosonic degree 
of freedom is retained as a representative of the continu- 
ous set of degrees of freedom. The function h(e) is cho- 



sen to be a constant in each intervall of the logarithmic 
discretization. The Hamiltonian is written in the new 
discrete basis and the resulting Hamiltonian is mapped 
onto a semi-infinite chain (Fig.[3J) with the electronic part 
H c \ coupling to the first site of the bosonic chain. Fi- 
nally the chain-Hamiltonian is numerically diagonalized 
via successively adding one site to the chain. The effec- 
tive Hamiltonian is treated on successive smaller energy 
scales by the renormalization group transformation 



H 



N+l — 



AH N 

A N+1 



£jv+i&J/+i&iV+i + tN (b^N+i + h.c. 



(26) 



The energies e„ and couplings between the elements of 
the chain are falling off as A~™. 

The bosonic NRG has been shown to give very accu- 
rate results for the spin-boson model |27ll28j . One of its 
strengths is the flexibility to handle a variety of models 
involving the coupling of a small subsystem to a bosonic 
bath. 

To detect possible phase transitions, we calculate the 
eigenvalue spectrum and the density-density correlation 
function (see below). In the limit of a — the exciton 
system and the bosonic degrees of freedom are completely 
decoupled. The coherent motion of the exciton is un- 
damped and we are in the delocalizcd phase. In contrast, 
in the case of J = the system is in the localized phase. 
The two phases (localized and delocalized) are connected 
by a quantum phase transition. Similar to the analysis in 
Ref. |2J| the phase diagrams of the exciton-boson models 
studied here can be obtained from the flow diagram of the 
lowest-lying many-particle levels. Another possibility is 
to calculate the density-density autocorrelation function 
C(uj) which shows a divergency at the phase transition. 

We calculate C[ui) for the dimer and trimer for differ- 
ent sets of parameters. This quantity is defined by 



1 

4-7T 



dte 



and probes the dynamics under equilibrium preparation. 
For the two-site model, C(u>) corresponds to the spin-spin 
correlation function of the equivalent spin-boson model 

1 r°° 

C(io) = —J dte Mt ([a z (t),a z (Q)}), 

with a z the z-component of the spin in this model. 

We calculate the density-density correlation function 
as sum of (5-functions in the Lehmann representation: 



C(«) = 5E 



UJ>0 



S (uj + e - e„) , 
(27) 
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The linear absorption and emission coefficient a(tu) for 
the donor site of the electronic system coupled to the 
bosonic bath under influence of an external laser field of 
frequency u) is given by Fermi's golden rule: 

a D {Lu) = 27r£|</ |ff£ rt | 0) | 2 5 (oj + E Q - E f ) , 
f 

(28) 

where Hp Crt is defined as 

H pert = C Dl C D0 +C ] m C Dl . (29) 

The term Hp eTt describes the excitation of an electron 
from the ground state Dq to the excited state D\ . It can 
be treated perturbatively as long as the probing photon 
energy is small. For the initial state we use the ground 
state |0) and |/) are all possible final states. 

The eigenenergies of H (Eq and Ef) and the matrix 
elements (0|-ff per t|/) are evaluated with the NRG for dif- 
ferent J and increasing coupling to the bosonic bath. To 
obtain a continuous curve for a (ui), the 5- functions ap- 
pearing in eq. EHI have to be broadened. Here we use 
the strategy discussed in Ref. that is replacing each 
(5-function by a gaussian on a logarithmic scale. On a lin- 
ear scale, this function is not symmetric around its center 
so that spectral weight in a D (uj) appears to be shifted to 
higher frequencies. 

RESULTS 
Phase Diagram 

Increasing J tends to delocalize the exciton. Since we 
have excluded the single-electron hopping in the Hamilto- 
nian, the dynamics can be restricted to the single-exciton 
subspace which maps onto a spin-boson model for two 
sites. For increasing coupling to the bosonic bath, the 
exciton localizes at a critical a c . We explore the phase 
diagram with localized and delocalizcd phases (connected 
by a Kosterlitz-Thoulcss transition) for the exciton-boson 
model with 2 to 6 molecular sites by calculating the crit- 
ical a c (J) as a function of J. Note that for a finite bias e 
(difference of the excitation energy on donor and accep- 
tor) there is no phase transition but a crossover from a 
delocalized to a more localized regime. 

The phase diagram for the spin-boson model was calcu- 
lated already in Ref. The critical a depends linearly 
on the matrix element J. It was noted that the exact 
value of a c has to be determined in the limit of A — > 1. 
We do not perform the extrapolation so that the critical 
a c is somewhat larger than the actual value. 

In Fig. we display the phase diagram of the multi-site 
exciton-boson model for 3,4,5 and 6 sites for both chain 
and ring geometry. 
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FIG. 6: Phase diagram of rings (solid line) and chains (dotted 
line) with 3,4,5, and 6 sites. The system is in the localized 
(delocalized) phase above (below) the phase boundary. 

The dashed and solid lines display the critical a c for 
the chain and ring, respectively. The critical coupling 
shows a linear behaviour similar as in the spin-boson 
model. For an even number of sites the ring has a larger 
critical a than the chain. For an odd number of sites, 
both curves cross at a certain value of J above which the 
opening of the ring will tends to delocalize the exciton. 

For a multi-site exciton-boson model with three or 
more sites, no quantum phase transition is observed as 
soon as the couplings Jy between neighbouring sites are 
different. To study the crossover from the delocalized to 
a more localized phase, we calculate equilibrium dynam- 
ical properties as discussed in the following. 

Equilibrium Dynamical Properties 

To study the dynamics of the electron transfer pro- 
cess in the one-exciton subspace we calculate the density- 
density correlation function C(uj) = ^ J_ °° e lwt C(i) dt 
with 

C(t) ~ (\cJ D a D {t) - a\aA{t),a\,a D {Q) - a A a A (0)] + ^ , 

(30) 

for the dimer (gn = 1/2, gA = — 1/2) and for the trimer 
(9b = 0,g A = -1,9d = !)• 

In the two-site case, the density-density correlation 
function is identical to the spin-spin correlation function. 
The correlation function shows the power-law behaviour 
for low frequencies up to w ~ T'. When a approaches 
a c , the slope in C(u) increases and the peak position (the 
characteristic energy and temperature scale T*) is shifted 
to lower energies, see Fig. [7| At the phase transition the 
correlation function is diverging. The correlation func- 



8 



100-, 




■ i ■ i 

0,00 0,02 0,04 

co/co 



c 



FIG. 7: Density-density correlation function C(u>) for the 
dirtier for different a and J = 0.4. For small frequencies, 
C(u>) is linear in u. The slope increases with increasing 
a — 0.1, 0.2, 0.3. For zero bias (e = 0) the phase transition to 
the localized phase is indicated by an infinite slope. 



Ji 




FIG. 8: Crossover temperature for a trimer ring with J2 = 
Jab = Jdb = 0.4 (tiangles), 0.5 (diamonds), 0.6 (circles) and 
various values of J\. The coupling a is set to 0.1 



tion shows an algebraic long time behaviour for T = 
and an exponential decay for finite T . 

In Fig.[3]we show results for the characeristic tempera- 
ture for a trimer ring by keeping Jab and Jdb constant 
and varying J\ = J da- {J da — corresponds to the 
chain.). The value of T* goes through a minimum at 
J\ = Jab = Jdb- The larger the difference the larger 
is the characteristic temperature. In Fig. [!|| we display 



1.0- 




FIG. 9: Crossover temperature for a trimer chain with J2 = 
0.5 and various values of Ji,. The coupling a is set to 0.1,0.2 
and 0.3 for the upper, middle and lower curve, respectively 

the characteristic temperature for an asymmetric three- 
site chain with J\ 7^ J 2 for increasing a. For large a 
and Ji = J 2 the characteristic temperature T* goes to 
zero indicating the phase transition. For the asymmetric 
chain no phase transition occurs and T* increases with 
the increasing difference of the matrix elements J\ and 

J 2- 

Absorption Spectrum 

For a — the spectrum of the two-site electron-boson 
model consists of four states. In the ground state, both 
electrons occupy the lowest level of donor and acceptor 
molecules, respectively. The system can be excited by a 
photon: D + A ~-» D* + A. If we now consider a finite 
J, the exciton is able to move to the acceptor and back 
(D* + D + A*). 

To calculate the absorption spectrum we choose the 
intitial state to be the ground state. The ground state 
is calculated with the NRG and depends on a and J. 
We set the energy difference between the ground state 
and the excited state to e = 0.75w c . For a = and 
J = the peak in the absorption spectrum is at u> — e = 
0.75. For increasing J peaks are at frequencies equal to 
the eigenenergies u> — e ± J. If now a is increased the 
two main peaks are broadened and shifted (see Fig. llOjl . 
The height of the peak at low frequencies increases with 
increasing a. 

The absorption spectrum for the trimer with a = 
shows peaks at e + 2 J and e — J for the ring geometry 
(with Jab — J ad — Jbd — J) and at e ± y/2J for 
the chain. The absorption spectra for various values of 
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FIG. 11: Absorption spectrum of a trimer ring for J = 0.2, 
€d = £A = 0.75u! c and a = 0.1, 0.2, 0.3 as a function of w. 



a are shown in Figs. UTI and for the ring and chain, 
respectively. 



CONCLUSION 

In this paper we studied the phase diagram, equi- 
librium dynamical properties and the linear absorption 
spectrum of Frenkel excitons in various ring and chain 
models with a coupling of the electronic degrees of free- 
dom to a bosonic bath. We used the numerical renormal- 
ization group method which allows to study the electron- 
boson and exciton-boson models in the full parameter 
regime. 



We studied in detail the phase diagrams of the multi- 
site electron-boson models in the subspace of one exci- 
ton. In the zero-bias case (all molecules have degener- 
ate HOMO and LUMO energies), increasing the value 
of a leads to a quantum phase transition between a de- 
localized and a localized phase. For the two-site case 
(dimer) the exciton-boson model can be mapped onto 
the spin-boson model for which the phase diagram is al- 
ready known. For more than two sites, the behaviour 
is more complicated and depends also on the geometry 
(chain vs. ring). 

The calculation of the density-density correlation func- 
tion allows to estimate the characteristic temperature for 
the crossover between delocalized and localized phase. 
This crossover temperature T* is zero for the localized 
phase and increases when the system goes to the delocal- 
ized phase. Increasing the difference between the various 
couplings of the model generally leads to a more delocal- 
ized behaviour. 

It would be interesting to compare our results to opti- 
cal experiments of small bio-engineered systems or quan- 
tum dots in which exciton transfer occurs. Further stud- 
ies are planned to evaluate the time dependent behaviour 
of excitons, to extend the system to larger rings and to 
include static disorder using time dependent and equilib- 
rium NRG to model systems like the LH II ring. 
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